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Abstract. The rehabihty of cluster lens reconstruction techniques based 
on weak lensing is studied in terms of the accuracy of their reproduc- 
tion of the total cluster mass as a function of distance from the cluster 
center. To do so, a variety of reconstruction algorithms is applied to 
synthetic lensing data created using a sample of 60 numerically mod- 
eled clusters, and the mass reconstruction is compared to the known 
deflector mass. The results can be summarized as follows: (1) Recon- 
struction algorithms which require integrations extending over the entire 
real plane yield unreliable results, because they give rise to boundary ef- 
fects which are hard to control; mass overestimates are more likely and 
more substantial in this case than underestimates. (2) Reconstruction 
techniques which avoid these boundary effects yield reliable lower bounds 
to the cluster mass. The tightness of such bounds depends on the size 
of the fleld, which can be extended synthetically to improve the results 
considerably. For the sample of numerical cluster models, the best lower 
bounds, achieved by combining synthetic fleld extension with non-linear, 
flnite-fleld reconstruction, decline from 100% to 80% of the true cluster 
mass going from the cluster center to an angular distance of 2.'5. The 
80% error bars of the lower mass bounds are ±10% to 15%. 
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1 Introduction 

It was realized by R. Webster as early as in 1985 that galaxy clusters could coherently 
distort images of background galaxies, and that information about the mass distribution 
in these clusters could be inferred from that distortion. After the detection of the faint 
blue galaxy population (Tyson 1988), flrst observational evidence for coherent distortions 
around galaxy clusters was found (Tyson, Valdes & Wenk 1990), and then attempts were 
made to systematically flt parametrized cluster mass distributions to the observed image 
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distortions (Kochanek 1990, Miralda-Escude 1991). Kaiser & Squires (1993) were the 
first to notice that, in the weak lensing regime, parameter-free cluster-lens inversions are 
possible on the basis of the observed image ellipticities. While Kaiser & Squires first 
suggested the application of their method to statistical cluster samples rather than to 
individual clusters except for very rich ones, several applications of their inversion technique 
to real data proved the feasibility of their method (Fahlman et al. 1994, Small et al. 1994, 
Kaiser et al. 1994b, c). 

While the overall agreement between the reconstructed mass maps and independent 
information derived from the galaxy- or the X-ray distribution demonstrated that the 
method worked successfully in principle, some intriguing puzzles were encountered. First, 
Fahlman et al. (1994) reported a total mass derived from lensing by the cluster MS 1224-|-20 
three times as large as the cluster's virial mass within an aperture of radius 2.'76 around 
the cluster center. This finding was qualitatively confirmed by Carlberg, Yee & EUingson 
(1994), who found the lensing mass 2.5 it 1.1 times larger than the virial mass. Fitting 
mass profiles to the weak shear observed in the field of CI 0024 -|- 1654, Bonnet, Mellier & 
Fort (1994) also found a surprisingly large mass for this cluster. Conversely, Squires (1994) 
reported that the cluster-inversion algorithm failed to detect any significant mass at all in 
the field of A 2163, which is particularly disturbing because this is the hottest and most 
luminous X-ray cluster in the sky. 

Recent theoretical studies aimed at removing sources of systematic error from the 
most straightforward cluster-reconstruction algorithm. As mentioned by Kaiser & Squires 
(1993) and discussed by Schneider & Seitz (1994), the surface-density map derived from 
coherent image distortions is invariant to a global linear transformation which expresses the 
fact that a matter sheet of constant surface density does not change the distortion pattern 
at all. This degeneracy can only be broken if independent information on the absolute 
surface mass density can be obtained, e.g., by observing the magnification of galaxy images 
behind clusters (cf. Broadhurst, Taylor & Peacock 1994, Narayan & Bartelmann 1994). 
The original method by Kaiser & Squires does only apply in the weak-lensing limit, i.e., 
if convergence and shear are small compared to unity. Two problems arise if this limit 
is not met, first, the observed image distortions then measure a combination of shear 
and convergence rather than the shear alone as in the case of weak lensing, second, the 
possible formation of critical lines by strong lenses causes a further degeneracy because 
then a given image distortion can arise outside and inside the critical curve, and they 
cannot be distinguished because the parity of images cannot be observed. Both these 
problems were solved by Seitz & Schneider (1994), who suggested an iterative algorithm 
capable of simultaneously reconstructing weak- and strong-lensing regimes (see also Kaiser 
1994). Another source of systematic error is caused by the fact that the most simple 
variant of the cluster-reconstruction algorithm requires a convolution to be performed on 
the entire real plane, while actual observational data are necessarily restricted to finite 
fields. If applied to finite fields, however, the cluster inversion yields a mass distribution 
with vanishing total mass in the field. A method to remove this systematic boundary effect 
was first suggested by Schneider (1994), and Kaiser et al. (1994c) discussed variants of 
this scheme. 

It is evident from the quoted references that this field of astrophysics is evolving at a 
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remarkable rate. This reflects the considerable interest attracted by the cluster-inversion 
technique, which is due to the facts that gravitational lensing provides essentially the only 
method to learn about the true mass distribution of clusters without referring to any 
equilibrium or stability assumptions, and that it was proved feasible to observe the weak 
lens effect by clusters with sufficient signal-to-noise ratio (see Kaiser, Squires & Broadhurst 
1994b for a description and discussion of methods of data analysis). 

This paper pursues the following question: Given "ideal" observational data, that is 
to say, observational data which have been cleaned of systematic observational errors due 
to the telescope or the CCD, and which are not contaminated by galaxies belonging to 
the cluster lens itself, what is the accuracy with which we can hope to reconstruct from 
weak lensing the total cluster mass as a function of distance from the cluster center? To 
answer this question, the weak lens effect of a sample of 60 numerically modeled galaxy 
clusters is simulated, and from the simulated background galaxy ellipticities the clusters 
are reconstructed using a variety of reconstruction techniques. To set the stage, I review in 
Sect. 2 the reconstruction algorithms whose application to the numerically modeled galaxy 
clusters is described in Sect. 3. In Sect. 4, I briefly discuss the influence of seeing and the 
unknown source redshift distribution, and the results are summarized and discussed in 
Sect. 5. 

2 Cluster mass estimators 
2.1 Basic relations 

This section reviews a number of different methods to estimate cluster masses from their 
weak lens effect. All these methods are physically based on the fact that the local properties 
of a lens, its convergence k, and shear 7, are related through the same scalar potential 



where a; is a two-dimensional position vector on the observer's sky. In the weak-lensing 
limit, the distortion of background galaxy images is determined by the shear, while the 
convergence is just the surface mass density in units of its critical value. Starting from the 
potential tp{x)^ the convergence is determined by the two-dimensional version of Poisson's 
equation, or 



where an index j preceded by a comma denotes partial differentiation with respect to xj. 
The components of the shear are likewise given by 




(2.1) 



^^i^) = = 2 t^'^^*^*^ + V',22(«)] , 



(2.2) 



7l(«) = 2 - V',22(«)] , J2{x) = tp,12{x) . 



(2.3) 



Equations (2.1) through (2.3) provide the complete basis for all cluster reconstruction 
techniques. For a detailed discussion of their derivation and their properties, see Schneider, 
Ehlers k Falco (1992, Chap. 5). 
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It will later prove useful to combine the shear components into the complex shear 7, 

7(a;) = 7i(a;) + i72(a;) . (2.4) 

Kaiser & Squires (1993) used the Fourier transform to show that k, and 7 are related by 
the convolution 

k{x) = -- [ d^x' Re [D{x - x')j*{x')] , (2.5) 



where the superscribed asterisk denotes complex conjugation, and the complex convolution 
kernel is given by 

V{x) = ^ ^ ^ . 2.6 

Equation (2.5) can be verified by inserting Eq. (2.1) into Eq. (2.3), obtaining a relation 
between 7 and and then inserting into that relation Eq. (2.5). 

2.2 Linear cluster reconstruction 

In the limit of weak lensing, I7I ^ 1 and «: ^ 1, the local ellipticity of galaxy images is 
a measure for the local shear 7(3;). The usual way to proceed is to measure the surface- 
brightness quadrupole Qjk{x) of each galaxy image, 

J d^x XjXkI{x) 

Qjk = — — — ' (2-7) 

where I{x) is the surface-brightness distribution of the image, and to form the complex 
(or two-component) ellipticity x, 

, . _ Qii — Q22 + 2i(5i2 /p, 

If the lens is weak, and if the lensed sources were intrinsically circular, the image ellipticity 
would be twice the shear, x = 27. In the most straightforward linear reconstruction 
algorithm, the ellipticity x is therefore taken as a local estimator of 27 and inserted into 
Eq. (2.5) to yield 

1 ^ 

where the area integral has been replaced by a sum extending over N galaxies, whose 
number density is n. 

In practice, a straightforward application of Eq. (2.9) is impaired by the fact that the 
lensed galaxies are intrinsically elliptical. The randomness of these ellipticities implies 
that the noise introduced by them has a white spectrum (see Kaiser & Squires 1993). The 
high-frequency contribution to the noise is therefore divergent, and this divergence must be 
suppressed by introducing a low-pass filter into Eq. (2.9). Intuitively, the divergent noise is 
caused by the singularity in the convolution kernel T>(x) for |a;| — )• 0, which adds divergent 
weight to the galaxy images very close to the point under consideration. Introducing the 
low-pass filter W(x; Ax) modifies Eq. (2.9) to read 
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1 ^ 

k{x) c:^-^J2w{x- x,-Ax) Re [V{x - x,)x*{xj)] . (2.10) 



27rn 



Since the only purpose of the fiher is to suppress the signal for any isotropic 

function W(x; Ax) is sufficient which tends to zero for |a;| — )• on a reasonable scale Ax^ 
and approaches unity for \x\ ^ Ax. An appropriate choice is (Seitz & Schneider 1994) 



2 \ / 2 

^ \ / X 



where Ax has to be chosen according to the spatial galaxy density, see below. Equation 
(2.10), supplemented by some choice of the smoothing kernel VF, provides a first simple 
estimator for the convergence k,(x) at a point x. 

2.3 Smoothing prior to reconstruction 

A reconstruction following Eq. (2.10) avoids the noise divergence caused by the intrinsic 
galaxy ellipticities, but still suffers from the randomness of the galaxy positions, which 
introduces "shot noise" into the reconstruction. A viable way to cure that is to smooth 
the observed data before applying the reconstruction algorithm to them, i.e., to construct 
a smooth ellipticity map from the data before inserting them into Eq. (2.5). For instance, 
such an ellipticity map can be determined at the points of a regular grid xjk by the 
weighted average 

X{Xjk) - — iv , (2.12) 

1^1=1 



where Gaussian weights wi appear appropriate, 

{Xjk - xi f 



wi = exp 



(2.13) 



Since Eq. (2.12) is a weighted average over the ellipticity data, it also reduces the noise 
introduced by the intrinsic galaxy ellipticities. Inserting x{x) into Eq. (2.5), we arrive at 
the alternative linear reconstruction equation 

zvr ^ — ' 

where a is the separation of grid points. 

An obvious choice for the smoothing length Ax, which occurs in the smoothing kernel 
W of Eq. (2.11) and in the Gaussian weights wi of Eq. (2.13), is the average distance 
between galaxies, 

Axq = ; (2.15) 



then, there will on average be one galaxy image within a smoothing disk. This is sufficient 
if the measured "signal" is strong, that is to say, if locally |x| — )• 1. It should suffice in 
that case to use the ellipticity of one or very few galaxies to estimate the local shear. If 
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the signal is "weak", the smoothing scale should be extended such that a smoothing disk 
contains a larger number of galaxies. However, what should be considered a weak signal 
depends on the intrinsic ellipticity distribution of the sources. To be specific, let us assume 
an intrinsic ellipticity distribution of the form 

Ps(|xs|)= T!'^~^''^^f^uK,^ ^ lxs|e[o,i], (2.16) 

TTxi [I - exp(-l/x^)J 

where Xr measures the scatter in the intrinsic ellipticity distribution. Clearly, the signal 
Ixl should be considered weak when |x| — )• Xr- Therefore, I adapt Ax to the strength of 
the signal such that 

. J 1 for Ixl - 1 . ^2_^7^ 



Axo \ 3 for Ixl Xr ' 

this means that a smoothing disk contains on average about ten galaxies when the signal 
is weak, and one galaxy when the signal is strong. A simple linear interpolation scheme 
which fulfils these conditions is 

Ax = Axo (3-2^^^) . (2.18) 

V 1 - Xr / 

To clarify the notation, 'linear reconstruction' without or with smoothing in the fol- 
lowing refers to the reconstruction equations (2.10) or (2.14), respectively, regardless of 
the fact that Eq. (2.10) already involves low-pass filtering the data. 

2.4 Iterative nonlinear reconstruction 

The linear reconstruction formulae (2.10) and (2.14) are valid in the limit of weak lensing. If 
this limit does not apply, several complications arise. First of all, Schneider & Seitz (1994) 
have emphasized that the only observable quantity accessible through the ellipticities of 
the galaxy images is the distortion (5, 

i = , ■ (2.19) 

(1 - Ky + 



which is related to the shear by 
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6' 



1±^1-\S\A ; (2.20) 



see also Miralda-Escude (1991). [The fact that only the combination S of k, and 7 is 
observable via ellipticities is most easily verified mapping a circular source with the inverse 
Jacobian of the lens mapping, A~ = diag[(l — k — j)~^ ,{1 — k-\-j)~^], and determining the 
ellipticity of the image.] The choice of the sign in Eq. (2.20) introduces a local degeneracy 
which is physically caused by the fact that a given distortion of the galaxy images can be 
produced on either side of the critical curve of a lens, because the parity of the images 
cannot be observed. The correct choice for the sign is the negative sign of the determinant 
of the local Jacobian matrix A. Two problems therefore arise immediately in the case of 
strong lensing, first, the distortion S has to be inferred from the ellipticity measurements. 
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second, the sign in Eq. (2.20) has to be chosen correctly. Both these problems have recently 
been solved by Seitz & Schneider (1994). 

They showed that the condition that the intrinsic galaxy ellipticities should be ran- 
domly oriented (i.e., x has random phase intrinsically) can be translated into the require- 
ment 

JV' r / •,2- 

(2.21) 



E 

J = l 



exp 



(X - Xjf 



1-Re [8{x)x] 







where the Gaussian factor was again introduced for smoothing the shot noise. This implicit 
equation for the distortion 8 can be solved iteratively, starting the iteration with (5*^"-' = x- 
The sum in Eq. (2.21) extends over all N' galaxies in a patch of the observed field large 
enough to contain a sufficient number of galaxies and small enough so that the distortion 
does not change significantly within the patch. Again, the size of the patch can be adapted 
to the strength of the signal, as described before. 

Having solved for the distortion field 8, Eq. (2.20) can be used to obtain an estimate 
for the shear. Since the convergence and the position of the critical curve are yet unknown, 
/€ = and det ^ > are assumed in a ffist step, yielding 



l-^l-\8{x)\^ 
8*{x) 



(2.22) 



The procedure now continues iteratively. In each step of the iteration, Eq. (2.5) is used to 
determine 



An) 



(x) = -- 



Re 



Vix - XjkYf 



— 1) 



{Xjk) 



(2.23) 



and from k*^"-* and the shear 7*^"-* 



^(-)(a,) = -^ 
vr 



jk 



V{x - Xjk)n'''^\xjk) 



the n-th iteration of the sign of the Jacobian can be found. 



Sl£ 



ign fdet^^")) = sign [(1 - k^'^^) 
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(n)\2 



(2.24) 



(2.25) 



The next iteration step starts again with the determination of 7 from Eq. (2.20), inserting 
the sign and n determined in the previous step. In practice, as long as the critical curve 
of the lens being reconstructed is not yet known, a further smoothing kernel should be 
introduced into Eq. (2.23) to smooth over the singularity which occurs in Eq. (2.20) if 
(5 — )• when the sign is negative. For an appropriate choice of this kernel, and for further 
details on this method, I refer to Seitz & Schneider (1994). 

The iteration is continued until a reasonable reproduction of the observed distortion 
field is achieved, which can be estimated by means of a x^ comparison between the observed 
and the reconstructed distortion field, given the intrinsic ellipticity distribution of the 
galaxy ellipticities. Again, Seitz & Schneider (1994) give a detailed description of the 
convergence properties of their iterative algorithm. 
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2.5 Boundary effects 

The three methods of cluster reconstruction described up to now aU involve convolution 
integrals extending over the entire real plane. Observational data, however, are limited 
to finite frames. If applied to finite data, the previously described techniques tend to 
construct spurious negative troughs in the surface-density maps. Several methods have 
been suggested so far to reduce or remove these effects. The one which is easiest to imple- 
ment was first mentioned by Seitz & Schneider (1994). They proposed to fit a spherically 
symmetric, analytic lens model to the distortion map obtained from the data, that is, the 
parameters of a lens model like 

--dei(x) = J^^^^f^ (2-26) 

are chosen such as to minimize the mean-square deviation between the distortion fields 
produced by the model and reconstructed from the data. The justification of this procedure 
is based on the expectation that the outer parts of clusters tend to become increasingly 
symmetric with increasing distance from the center. Therefore, the adaptation of the lens 
model is performed using only the fraction of the field outside a specified radius Xflt, which 
can be on the order of arc minutes. This procedure can be used to artificially extend 
the field into regions where k, and 7 are well below any measurable value. I refer to this 
technique in the following as "field extension" . 

Other more involved methods to remove the boundary effects introduced by the finite 
size of the field start with Kaiser's (1994) realization that the gradient of k, can be written 
as a linear combination of derivatives of 7, 

Vk=( ^1-1+ ^2,2 ^ (2.27) 
V72,i - 71,2 / 

so that the difference in k, between any two points, Xq and a;, can be found by integrating 
the components of V/€ along arbitrary lines connecting the two points, 

j^X 

k{x) = k{xo)+ dl-VK. (2.28) 

J Xo 

To obtain an estimate of n{x) at a given point a;, one can choose Xq on the boundary of 
the field, integrate towards a;, and then average the result over all boundary points. If the 
boundary is sufficiently distant from the center of the cluster, h{xq) can safely be assumed 
to vanish; otherwise, the procedure estimates the difference 

k(x) - j) d(f)K[xo((f))] , (2.29) 

where Xo((f>) parametrizes the boundary of the field. A variant of this technique to re- 
construct K, was suggested by Schneider (1994). Another way to use Eq. (2.27), which is 
numerically easier to implement, was suggested by Kaiser et al. (1994c); they proposed 
to average the starting points Xq of the line integrals over the entire observed field. As 
shown below, this essentially leads to modification of the convolution kernel of the previous 
reconstruction formulae. 
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Suppose the boundary of the field be given by a closed curve c(^), parametrized by 
the continuous parameter (f> G [0,27r]; further, let R{(f>) = |c(^)| be the distance of the 
boundary from the center of the field, and A be the area enclosed by c. Without loss of 
generality, let the target point, i.e. the point where k, should be reconstructed, be at the 
coordinate origin. Then, transforming coordinates {xi^X2) to (t,^), where 



xi = tR cos • 



X2 = tR sin ( 



(2.30) 



radial lines connecting the boundary with the origin are given by l{t,cf)) = t(cos^,sin^) 
with t G [0, 1], and the average of Eq. (2.28) over A reads 

/c(0) - R = 

- — / R {ct>)dct> I tdt I dt' [cos ^(71,1 + 72,2) + sin ^(72,1 - 71,2)] • 
^ Jo Jo Jo 

K, denotes the average value of n within A. Since the derivatives of 7^ have to be taken 
along the paths, they depend on t', and therefore the t-integration can be performed. 
Renaming t' to t afterwards, we obtain 



2A 



27r 



R^(f))d(f) / dt(l-f) [cos (^(71,1 +72,2) + sin (^(72,1 -71,2)] . 



(2.32) 



Transforming back to cartesian coordinates (2.30), and inserting cos ^ = xi/x, sin^ = 
X2/x^ Eq. (2.32) can be written 

k{0) - K = ^ J dxi J dx2 ^1 - (a;i7i,i + ^272,1 + a;i72, 2 - 2:271,2) • (2.33) 



Using the relations 



1 



R' 



2R 



Rxi + Rx-^ 



d 
dxi 
d 

dx2 \ x^ I 

where R denotes {dR/ dcf)), Eq. (2.33) can be integrated by parts to yield 
/c(0) - R = 



1 = — RX2 - Rxi 



(2.34) 



1 



d^x 



R^ RR 

— (71 cos 2(f) + 72 sin 2(f)) + —j- (71 sin 2(f) - 72 cos 2(, 



X 



(2.35) 



Shifting the target point into an arbitrary position x within the boundary c(^), we finally 
obtain 



K,(x) — H 



J J d^x' i^^R^Re[V(x-x')-f*(x')]-RRlin[V*(x-x')-f(x')]^ . 



(2.36) 
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Apart from changes of notation, Eq. (2.36) agrees with Eq. (7) of Kaiser et ah (1994c). Note 
that (2.36) reproduces (2.5) when the boundary is extended to infinity since then {R^ / A) — > 
(l/vr) and i? — )• 0. Ahke the previous area integrals, Eq. (2.36) can be approximated by 
a finite sum over galaxy images, but now it extends over a finite area only, thus avoiding 
boundary effects. Note, however, that R is the distance between the target point x and 
the boundary via the point a;', and R is the derivative of this distance with respect to the 
polar angle around the target point. If x is close to the boundary, R can become very 
large. In addition, if a galaxy happens to be very close to the target point, a small change 
of the galaxy position can substantially change R and R. Both these effects can introduce 
loud noise due to the random positions of the galaxies, if Eq. (2.36) is approximated by a 
discrete sum over individual galaxy images. It is therefore essential to smooth the data, 
e.g., according to Eq. (2.12), before inserting them into Eq. (2.36). Even after smoothing, 
it is not possible to reconstruct n with Eq. (2.36) out to the field boundary; however, 
simulations show that Eq. (2.36) provides a successful reconstruction mechanism even close 
to the boundary. 

Remember that the purpose of this study is to assess the reliability of the total recon- 
structed cluster mass. Although the finite-field kernel removes the systematic boundary 
effect from the reconstruction, it yields k, — R rather than k. If we integrated k, — R out to 
the field boundary to determine the total enclosed mass, the result would vanish again, by 
definition of and the only advantage of the finite-field reconstruction would be that the 
total mass tends to zero not because of systematic errors, but because the underlying R of 
the cluster is unknown. For a mass estimate based on the finite-field reconstruction to be 
reliable, it is therefore necessary that the mass determination be limited to an area smaller, 
and preferably much smaller, than the observed field, in order to make R much smaller 
than K. In any case, a reconstruction according to Eq. (2.36) provides a lower bound to the 
total cluster mass, and this bound will become tighter when the observed field is enlarged 
- see Sect. 3.3 below. 

The modified convolution kernel in Eq. (2.36), which I will call "finite-field kernel" in 
the following, can also be inserted into Eq. (2.23) to replace the kernel previously used in 
the non-linear reconstruction algorithm. Then, Eq. (2.23) reconstructs the n-th iteration 
of K, — R instead of k, itself, however avoiding systematic boundary effects. If inserted 
into Eqs. (2.24) and (2.25) to iterate the critical curve and the sign of the Jacobian, this 
introduces an inaccuracy. In addition, the kernel used in Eq. (2.24) to iterate the shear is 
still the kernel which applies to the entire real plane instead of the finite field. However, 
both these inaccuracies are negligible. If we are close to a critical curve, «: ^ 1, hence R 
is much smaller than k, in this case if the field is not too small. Furthermore, the only 
purpose of Eq. (2.34) is to iterate the shear for the determination of sign(det^). Since 
this sign will only change from -|-1 to —1 very close to the cluster center, boundary effects 
introduced by the truncation of the convolution do not matter for that purpose. I will 
apply this finite-field nonlinear reconstruction to model clusters in the following section. 

2.6 Global statistical measures 

Instead of reconstructing the entire two-dimensional surface-density map of the cluster, 
the more modest goal of this paper is to estimate masses from weak lensing. A simple way 
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to determine the average convergence enclosed by a boundary starts by using Poisson's 
equation to write n as the divergence of a vector ol = V^/', 

K = ^V-cx- (2.37) 

OL is the deflection angle caused by the lens. With the two-dimensional version of Gauss' 
theorem, the area integral over n within a boundary curve c{cf)) can be written 

(fxK{x) = -(td(i)a-n, (2.38) 
A 2 / 

where n is the outward-directed vector orthogonal to the boundary's tangent direction, 
n = (c2, ~ci), and thus 

k{x) = - (t d(f) (aiC2 — a2Ci) . (2.39) 
A 2 / 

Integrating by parts, Eq. (2.39) can be transformed to 

J Sx k(x) = - j) d(f) (a2Ci — aiC2) . (2.40) 

If we now specialize c(^) to a circle, i.e., ci = i?cos^, C2 = -Rsin^, we can write 

ai = ciai^2 — C2ai,i , ^2 = Cia2,2 — C2a2,i • (2-41) 
Moreover, the deflection angle a is the gradient of hence 

aj,k = 4^,jk = 4^,kj , (2.42) 

and, using Eqs. (2.2) and (2.3), 

+ 71 5 "2,2 = « -71 5 "1,2 = 72 • (2.43) 
Inserting this together with (2.41) into (2.40), we arrive at the result 

^=jj d^xK = («:) + (7t) , (2.44) 

where 7t is the tangential component of the shear along the boundary, 

7t = 71 cos 2^ 72 sin2^ , (2.45) 

and the angular brackets in Eq. (2.43) denote the average over the boundary curve. 

Since (k.) > along any boundary, Eq. (2.44) provides a lower limit to the mass 
enclosed by the boundary, 

Hr) > (7t) • (2.46) 

On the other hand, since 



2 r , , , 

K{r) = ~ r'dr' (/€)(r') 
f Jo 



(2.47) 
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(2.48) 



Eq. (2.44) can be transformed to 

2(7t) ^ dR{r) 
r dr 

Integrating Eq. (2.48) between radii ri and r2 > ri, we finally find 

C{r,,r2) = R{n)-R{n,r2) = -^^ / — (7t) , (2.49) 

and the integral can be written as an area integral, 

C(ri,r2) = -- / d^xRe[V{xyr{x)] . (2.50) 
vr r| - rf J 

The function C{ri^r2) has been derived by Kaiser et al. (1994a), who called it the (- 
statistic. Since «:(ri,r2), the average convergence in an annulus bounded by radii ri and 
r2, is non-negative, C{ri^r2) provides a lower bound to the average convergence within ri. 
In principle, Eq. (2.46) could already be used to provide a lower bound on the average 
convergence within a circular aperture, but it is observationally easier to approximate the 
area integral in Eq. (2.50) by a sum over all galaxies within a suitably sized annulus and 
use Kaiser's C-statistic lower bound on «:(ri). 

If applied to a set of nested apertures. Kaiser's C-statistic can be used to estimate lower 
bounds on the masses enclosed by the apertures. Consider a set of radii r^, 1 < z < n. If 
( has been measured in all rings with ri < rj^ Eq. (2.49) can be written 

f^i f^ij — Cij 7 (2.51) 

where Ri = R{ri) and Rij = R{ri,rj). The equations (2.51) are not independent since they 
describe in essence the cluster masses in rings, which must add up to the same mass inside 
a given radius irrespective of how the radius is decomposed into rings. The masses Mij in 
rings ri < r < rj are 

M,, = A,, R,j , = {r) -rjy , (2.52) 

and they have to satisfy the relations 

= Mk-i,k , M, = Y, Mk-i,k , (2.53) 

k=i+l k=l 

where Mo,i = Mi. For convenience of notation, we write for the masses inside adjacent 
rings, nik = M^-i^k- Then, inserting (2.52) and (2.53) into (2.51), we get 



rrik- Y rrik = A^jQj ; (2.54) 

k=l k=i+l 



a set of (n/2)(n — 1) equations for the n unknowns nik- Of these equations, however, all 
equations for fixed j are redundant, and thus (2.54) provides essentially a set of n — 1 
equations for n unknowns. The fact that there is one equation less than unknowns reflects 
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the scaling invariance of n. A reasonable choice is to set m„ = 0, i.e., to assume that there 
is no mass in the outermost aperture. Then, Eqs. (2.54) yield for mi, 

mi = Ai ^^,nCi^n-A, ,Ci,n-i ^2_gg^ 
Keeping i = 1 fixed, Eqs. (2.54) simplify to 



rrik = —T'rni - Ai^jd^j , (2.56) 

k=2 



and the latter set of equations can then be used to determine the masses Mj enclosed by 
the radii rj, 

= + , (2.57) 

for J > 2. The masses Mj are still a lower bound to the actual masses in circles around 
the cluster center, because m„ = was chosen, but this lower bound will be tighter than 
( alone. 

Note that the restriction to circular boundaries made before Eq. (2.41) simplifies the 
equations, but is by no means necessary; any shape of a closed boundary would be appli- 
cable. In practice, this means that the boundaries chosen could be adapted to the shape 
of the cluster under consideration, i.e., the boundaries could be shaped like the X-ray 
isophotes or the galactic isopleths of the cluster. 



3 Application to cluster models 



3. 1 Sample of model clusters 

The numerical cluster models that will be used in the following are taken from the sample 
described and used in a sequence of earlier papers (Bartelmann & Weiss 1994, Bartelmann, 
Steinmetz & Weiss 1994, Bartelmann 1994). Briefly, they were produced by iV-body 
simulations starting from CDM initial conditions which were normalized to the COBE 
quadrupole mesurement of the CBR, for the cosmological parameters i?o = 1, A = 0, and 
Hq = 50 km/s/Mpc. In total, 13 clusters were simulated, and their three-dimensional 
particle distribution was stored at about ten time steps per cluster between redshifts 1 and 
0. The projections of these three-dimensional models along the three independent spatial 
directions can serve as independent cluster models for the purposes of this paper. 

From this numerical cluster sample, I select 60 cluster models with redshifts in the 
range Zd G [0.1,0.4] with a large spread of their maximum convergence Kmax £ [0.4,1.2] 
with a median of Rmux = 0.8. The distribution of Kmax is displayed as a histogram in 
Fig.l. 

The clusters are highly irregular in shape, and many of them show substructure and 
mergers of smaller mass condensations. The flelds containing the clusters are square-shaped 
and have a uniform angular side length of 5'. 
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Figure 1.- Histogram of the maximum conver- 
gence Kmax of the 60 cluster models used in 
this paper. The median of the distribution is 

^ma,x — 0.8 



3.2 Determination of imaging properties 

To simulate the distortion of background galaxy images caused by the cluster models, 
the cluster fields were populated with randomly placed galaxies such that their averaged 
number density was n = 35 (arc min)~"^, corresponding to an average of N = 875 galaxy 
images per cluster field. Their intrinsic ellipticities Xs were drawn from a distribution of 
the form (2.16), with width Xr = 0.3; for this choice, see Miralda-Escude (1991) and Tyson 
& Seitzer (1988). With the known convergence and shear fields of the cluster models, the 
intrinsic ellipticities can be transformed to image ellipticities with 



X 



Xs+2g + g\* 
l + \g\^ + 2Re{gxt) 



where g is the combination 



9 



7 



1 



(3.1) 



(3.2) 



Once the image ellipticities are found, they can be inserted into the reconstruction equa- 
tions derived in the previous section. 

For each method applied to the cluster sample below, I determine the average recon- 
structed convergence within circles centered on the cluster center and normalize it by the 
true average convergence of the input model; the results therefore give the fraction of the 
reconstructed mass in units of the total mass as a function of distance from the cluster 
center. The error bars in the following figures show the range containing 80% or 48 of the 
cluster models, and the curves show the average. 

In the following, I apply all the reconstruction methods described above to the numer- 
ical cluster sample, viz. the linear reconstruction algorithm without and with smoothing of 
the input data, the nonlinear iterative reconstruction without and with field extension, the 
finite-field convolution kernel, either inserted into the linear or the nonlinear reconstruc- 
tion formulae, without and with field extension, and finally the global mass reconstruction 
within nested circular apertures, again without and with field extension. Where applicable, 
the data are smoothed with the adaptive smoothing length given in Eq. (2.18). If the field 
is synthetically extended, it is extended to twice its side length, i.e., to 10'. The number 
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of grid cells in the case of the two-dimensional reconstruction techniques is usually 32^, 
which is increased to 64^ when the field is extended in order to keep the spatial resolution 
the same. All the results given below concern the reconstructed mass fraction, i.e., the 
ratio between the reconstructed mass within a certain radius Mj-ecir) and the true mass 
^true(?") within the same radius. For convenience of notation, this ratio is abbreviated by 

M,ec(r) 

We can roughly estimate the expected behaviour of /(r). For all reconstruction methods 
involving infinite integrations, /(r) will go to zero towards the field boundary. For the 
finite-field reconstructions, assume for simplicity that the cluster is approximately isother- 
mal at sufficient distance from a possible core, hence Ktrue (V'") ^true Since 
n — K, is taken as an estimate for n in this case, /(r) will behave like 

— 2 

/(r) ~ 1 - ^ ~ 1 - Cr , (3.4) 
where C oc i.e., /(r) will decline linearly from ~ 1 at small r outward. 



3.3 Mass reconstruction from two-dimensional density maps 

As a first example. Fig. 2 displays the result of applying the unsmoothed linear recon- 
struction formula (2.10) to the simulated image data produced by the numerical cluster 
sample. 
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Figure 2.- Linear reconstruction according to 
Eq. (2.10), without smoothing the input el- 
lipticity data. Close to the cluster center, 
the mass is reasonably reproduced, but the 
boundary effects of this method require the 
total mass within the observed field to vanish, 
hence the reconstructed mass fraction tends 
to zero towards the field boundary 



The two hatched regions in the figure shows the 80% error interval for the mass recon- 
struction, together with the average over the sample. It is apparent from the figure that 
even the most straightforward application of the linear cluster reconstruction algorithm 
produces reasonable results close to the cluster center; at r = 0.25', f(r) = 0.90 it 0.2. 
With increasing radius, however, /(r) tends to zero. The reason for this behaviour is that 
the applied reconstruction technique approximates a convolution which should extend over 
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the entire real plane by an area integral over a finite field, and therefore the method re- 
constructs a mass distribution with vanishing total mass within the field. All the more 
elaborate versions of the reconstruction formula show this behaviour as long as they involve 
an integration over the entire plane. Panels (a), (b) and (c) of Fig. 3 show examples. 

Panel (a) is a reproduction of Fig. 2, added for comparison. Panel (b) shows the result 
obtained with Eq. (2.14) after smoothing the ellipticity field according to Eq. (2.12). The 
smoothing of the input data reduces the signal and therefore leads to a lower average /(r) 
for r — )• 0, but the scatter in /(r) is reduced; again taken at r = 0.25', f(r) = 0.73 it 0.13. 
The nonlinear reconstruction following Eqs. (2.19) through (2.25), the results of which are 
displayed in panel (c), exhibits the same qualitative behaviour, slightly further reducing 
the scatter; /(r) = 0.75 ± 0.09 at r = 0.25'. 
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Figures.- Comparison of several variants of the reconstruction algorithm. Panel (a): reproduction of 
Fig. 2, i.e., linear reconstruction without smoothing the ellipticity data; panel (b): linear reconstruction 
after smoothing the ellipticity data; panel (c): nonlinear reconstruction; panel (d): reconstruction with 
the finite-field kernel of Eq. (2.36). From (a) to (b), part of the signal is lost due to smoothing, but the 
errors towards the center decrease. The nonlinear reconstruction (c) further reduces the error close to the 
center. Panels (a) through (c) clearly show that the boundary effects require a vanishing total mass in the 
field. The finite-field reconstruction (d) removes this boundary effect, but increases the error 



Besides the unwanted feature of reconstructing a vanishing total mass, these results 
show that the error in the straightforward linear reconstruction formula (2.10) can be re- 
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duced by about a third when the input data field is smoothed prior to the convolution, 
and by about a half when the additional modifications required by the nonlinear recon- 
struction algorithm are included. The amount of the latter reduction, which comes from 
a more reliable reconstruction of the clusters with high Kmax in our numerical sample, 
depends of course on the distribution of peak values for n in the actual cluster sample. For 
comparison, panel (d) of Fig. 3 shows the result of the finite-field reconstruction following 
Eq. (2.36), i.e., a reconstruction which does not involve an infinite integration. As was to 
be expected, the mass-reconstruction function /(r) does no longer tend to zero when r 
approaches the field boundary, but remains finite. 

One could attempt to employ the scaling invariance of the distortion field b to rescale 
the convergence field by the global transformation 



(1-A) + Aa 



(3.5) 



in order to correct for the systematic boundary effect which forces the total mass in the 
field to vanish. Since the decay of /(r) for r going towards the field boundary arises in 
part because the simpler variants of the reconstruction method create spurious negative 
troughs in the convergence field, a plausible rescaling of n would be to choose A such that 
/€ is non-negative everywhere. As an example, we display in Fig. 4 the results of such a 
rescaling of the convergence maps obtained from the linear reconstruction equation (2.10). 
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Figure 4.- Linear reconstruction without 
smoothing the ellipticity data, scaled with the 
global transformation (3.2) such that «; > 
everywhere in the field. Panel (a): entire 
sample; panel (b): Kmax < 0.8, panel (c): 
'tmax > 0.8. While the rescaling yields rea- 
sonable results for the high-^max subsample, 
it tends to substantially overestimate the to- 
tal mass for large radii 



Panel (a) of the figure shows the results for the entire cluster sample, panel (b) for 
the low-/€ clusters (/Cmax ^ 0.8), and panel (c) for the high-/€ clusters. While the results for 
the high-/€ clusters are fairly accurate, /(r) G [0.87. . . 1.28] at r = 2.5', the errors for the 
low-/€ clusters become substantial at this radius, /(r) = [0.94. . . 1.75]. The reason is that 
this rescaling procedure exaggerates the necessary correction. If the signal, i.e., the shear- 
induced ellipticities of the images becomes small compared to the intrinsic width of the 
source ellipticity distribution, substantial fiuctuations in the convergence map arise which 
are not systematic like the spurious negative troughs, but random. Since this noise has a 
white spectrum, positive and negative fiuctuations of that kind should cancel on average, 
but the rescaling described here lifts the convergence map so far that «: > everywhere. 
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Therefore, positive fluctuations in k, can no longer cancel out, and the rescaled convergence 
map overestimates the total mass. This eff'ect is strongest for low-n clusters since there 
the overall signal is weaker compared to the intrinsic ellipticity. 

As mentioned before, the convolution kernel of Eq. (2.36), which was derived to avoid 
systematic boundary effects, can straightforwardly be inserted into the nonlinear recon- 
struction algorithm, which should result in a better mass reconstruction for small radii. 
Panel (a) of Fig. 5 shows the results obtained with that technique. 
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Figures.- Nonlinear reconstruction using the 
finite-field convolution kernel of Eq. (2.36). 
Panel (a): using only data in the original field; 
panel (b): after extending the field synthet- 
ically. The nonlinear reconstruction reduces 
the error close to the cluster center (cf. panel 
(d) of Fig. 3), and the field extension improves 
the reconstruction result considerably 



As can be seen by comparing panel (a) of Fig. 5 to panel (d) of Fig. 3, the nonlinear 
reconstruction algorithm reduces the error for small r, while leaving the results for larger 
r basically unchanged. Taken again at r = 0.25', f{r) G [0.59... 1.1] for the linear and 
/(r) G [0.64 . . . 0.94] for the nonlinear reconstruction. 

Although Eq. (2.36) avoids the unwanted negative troughs, it has the disadvantage 
that it does not reconstruct the convergence fleld, but rather n{x) — i.e., the convergence 
reduced by n averaged within the fleld. (Note that this is not a transformation of the type 
(3.2).) It therefore seems promising to reduce «: by extending the fleld. In panel (b) of 
Fig. 5, I show the results for /(r) obtained in the following way. First, the data fleld is 
extended according to the prescription (2.26), i.e., by adapting a spherically symmetric, 
parametrized lens model to the distortion data and using this lens model then for creating 
synthetic distortion data outside the 'observed' fleld. Then, the nonlinear reconstruction 
technique with the flnite-fleld convolution kernel of Eq. (2.36) is applied to the extended b 
fleld. Especially in comparison with panel (a) of Fig. 5, panel (b) shows that this procedure 
is successful; it not only reduces the noise of the mass reconstruction but also corrects for 
part of the decline for increasing r apparent in panel (a). 

3.4 Lower bounds to the radial mass profile 

For the sole purpose of reconstructing masses within apertures, the construction of two- 
dimensional mass maps is not necessary. I therefore also apply Eq. (2.57) to the synthetic 
image data; the results are displayed in Fig. 6. 

The comparably small error bars apparent in Fig. 6 are due to the fact that, in de- 
termining the (ij required in Eq. (2.57), the image data are averaged within much larger 
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Figure 6.- Lower bounds to the mass enclosed 
in circular apertures. Panel (a): using only 
the observed field; panel (b): after syntheti- 
cally extending the field. Especially after field 
extension, the mass reconstruction close to 
the cluster center is very good, but it drops to 
~ 60% of the mass close to the field boundary 

parts of the observed field. Since the boundary condition m„ = entered the derivation 
of Eq. (2.57), the result is a lower bound to the actual mass. It underestimates the true 
mass by more than 50% close to the boundary of the field; see panel (a). Again, the re- 
sults improve if the data field is extended with the help of a parametrized, symmetric lens 
model. Panel (b) of Fig. (6) shows the result for the lower bound to the radial mass profile 
according to Eq. (2.57), obtained after extending the distortion field. The field extension 
increases the recontructed mass fraction from 40% to 60% close to the field boundary, 
and it also improves the accuracy of the reconstruction towards the center of the field. 
At r = 0.25', /(r) = 0.95 ± 0.15 for the original data field, and /(r) = 1.0 ± 0.1 for the 
extended data field. 



4 Influence of seeing and tlie source redsliift distribution 

Seeing reduces the signal because it tends to circularize images. The effect of seeing can be 
approximately quantified in the following way. Suppose the surface-brightness distribution 
of the source is Gaussian, 

/(a.)ocexp(^--i-^j , (4.1) 

where a and h are the semi-major and semi-minor axes of the image ellipse. The effect of 
seeing is modeled by convolving I{x) with a Gaussian of width s. It is then straightfor- 
ward to show, most easily with the convolution theorem for Fourier transforms, that the 
ellipticity \x\' of the image affected by seeing is related to |x| by 





X 




1 + ^2(1 




xl) 



(4.2) 



where a is the width of the seeing disk in units of the intrinsic source width 6, a = (s/b). 
Since the reconstructed convergence is linear in x within the largest fraction of the observed 
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Figure?.- Influence of seeing on the elliptic- 
ity of galaxy images; their surface-brightness 
proflle was assumed Gaussian. The flgure 
shows the ellipticity after seeing, divided by 
the true ellipticity, as a function of the width 
of the seeing disk in units of source widths, 
cr = (s/b). The effect of seeing becomes 
stronger for weaker signal, as expected. Sig- 
nals of moderate ampitude, |x| — 0.5 say, are 
reduced by ~ 50% if cr ~ 1.4 



cluster field, Eq. (4.2) gives an estimate for the effect of seeing on the surface- density maps. 
Fig. 7 displays (IxI'/lxl) dependence of a for four different values of 

Clearly, a very strong signal, |x| — )• 1, is hardly affected by seeing at all; it requires a 
width of the seeing disk of 



a 



Ixl 



(4.3) 



in order to reduce the signal by 50%. A signal of moderate strength, say |x| = 0.5 . . . 0.6, 
is reduced to one half if the seeing is roughly 1.5 times the source size. In principle, it is 
possible to correct for the decrease of the signal due to seeing using Eq. (4.2). 

Another systematic source of error is the unknown redshift distribution of the sources. 
Assuming a certain source redshift zl. different from their true redshift Zg leads to an 
overestimate of the cluster mass if Zg > zl.^ and conversely. This can be quantified by the 
geometrical factor entering the conversion between surface mass density and convergence. 



K(zd,Zs) OC 



(4.4) 



where Dd,s,ds are angular-diameter distances from the observer to the cluster and to the 
source, and from the cluster to the source, respectively. Fig. 8 displays the function 



r(zd;zs,z'j 



K{zd,Zs) 



(4.5) 



i.e., the ratio between the convergences for the true and the hypothesized source redshift, 
respectively. The definition implies that r is the factor by which the cluster mass is 
overestimated if z'^ is assumed in place of the true source redshift. For the figure, 2;^ = 1 is 
kept fixed, and curves for four cluster redshifts Zd G {0.1,0.2,0.3,0.4} are shown. 

The figure shows that the uncertainty due to the unknown source redshifts is negligible 
if the cluster is at a redshift significantly lower than the sources; for Zd = 0.2, r G [0.89, 1.06] 
for Zs G [0.6,2.0]. For high-redshift clusters the error can be substantial; r G [0.33,1.34] 
for Zd = 0.5. For moderate-redshift clusters, however, the uncertainty in the source red- 
shifts translates to an uncertainty in the cluster mass on the order of 10% . . . 20%. If the 
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Figure 8.- Uncertainty caused by the unknown 
source redshift distribution, quantified by the 
function r defined in Eq. (4.5); r is the factor 
by which the cluster mass is overestimated if 
the sources are assumed at redshift = 1 
while their true redshift is z^. The result de- 
pends on the cluster redshift as indicated; 
for clusters at low redshifts, the uncertainty is 
negligible, and it increases for increasing z^ . r 
is very insensitive to the cosmological param- 
eters J7o and /Iq ; for the plot, an Einstein-de 
Sitter model was assumed 

sources are distributed in redshift, r has to be integrated over the source-redshift interval, 
weighted by the redshift distribution function. The curves in Fig. 8 were calculated for the 
parameters of an Einstein-de Sitter universe, but they are very insensitive to changes in 
the cosmological model. 

5 Summary 

A number of different cluster reconstruction techniques based on weak lensing were re- 
viewed and applied to a sample of numerically simulated cluster mass distributions. The 
following methods were discussed: 

1. Linear reconstruction without smoothing the observed galaxy ellipticities (Kaiser & 
Squires 1993). This method is designed for weak lens effects, i.e., for n ^1 and 7 ^ 1. 
It involves low-pass filtering of the signal to reduce the otherwise infinite noise due to 
the random ellipticities and positions of the galaxies. 

2. Linear reconstruction after smoothing the galaxy ellipticities (Seitz & Schneider 1994). 
This method is a variant of (1.) and still operates in the weak-lensing regime. 

3. Non-linear iterative reconstruction (Seitz & Schneider 1994). This technique proceeds 
beyond the weak-lensing limit, and therefore has to control the critical curves of the 
cluster. 

4. Finite-field reconstruction (e.g., Schneider 1994, Kaiser et al. 1994c), which replaces 
the convolution kernel operating on IR required by methods (1.) through (3.) to 
convert shear into convergence by a kernel which operates on a finite field only. This 
replacement avoids boundary effects which reduce the total mass within the observed 
field to zero. This finite-field kernel can be used by methods (2.) and (3.), regardless 
of non-linearities in the lens effect. 

5. Synthetic field extension (Seitz & Schneider 1994), which was suggested as another 
means to reduce the boundary effects of methods (1.) through (3.). It proceeds by 
adapting a spherically symmetric lens model, parametrized by central convergence, 
core radius, and convergence-profile slope, to the outer parts of the observed ellipticity 
data, and then synthetically extending the measurements beyond the field boundaries. 
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Again, this field extension is independent of the reconstruction technique chosen and 
can be combined with either. 
6. Mass reconstruction within circular apertures (based on the C-statistic by Kaiser et 
al. 1994a). This method employs the fact that, by Gauss' theorem in the plane, 
the average convergence within a boundary can be expressed by a line integral over 
observable quantities along the boundary. 
The results described above can be summarized in the following way. 

a. The accuracy of the mass reconstruction close to the cluster center is acceptable even 
for the most straightforward application of the linear reconstruction algorithm. The 
error in the function /(r) defined in Eq. (3.3) is ±0.2% for method (1.), and can be 
reduced by methods (2.) and (3.) to ±0.13 and ±0.09, respectively. The improvement 
achieved by method (2.) is a result of smoothing the observed data prior to use, and 
the further improvement by method (3.) is due to the more accurate reconstruction 
of the centers of clusters with central convergence Kmax — 1- 

b. A substantial disadvantage of methods (1.) through (3.) is that they approximate an 
integral which should cover IR by an integral over the observed field only. As a result, 
they reconstruct a mass distribution in the field with vanishing total mass by producing 
spurious negative troughs in the convergence map close to the field boundary. A 
possible cure would be to employ the global scaling invariance (3.5) to scale k, such 
that it is non-negative everywhere in the field. This, however, removes some negative 
fluctuations in k, arising from the statistical fluctuations in the intrinsic ellipticities 
of the sources which should in the ideal case cancel against positive fluctuations. 
Correcting such that k, > everywhere in the fleld thus overestimates the cluster 
mass especially if the cluster's lens effect is weak, because if k, is small, the statistical 
fluctuations are relatively stronger compared to the true k, than if the lens effect is 
stronger. While the rescaling to non-negative k, works quite well for clusters with high 
central convergence, Kmax ^ 0.8 say, it can produce a substantial overestimate of the 
cluster mass, by up to a factor of ~ 1.8 in our simulations, depending on the size of 
the fleld. 

c. The method to correct for this boundary effect employed here, based on Eq. (2.36), is 
successful, but it suffers from a number of different shortcomings. First, it does not 
reconstruct k, directly, but rather the difference k, — i.e., the convergence reduced 
by the average convergence in the fleld. However, R can be reduced by enlarging 
the fleld. Second, the flnite-fleld convolution kernel of Eq. (2.36) is susceptible to 
noise especially close to the fleld boundary and if unsmoothed observational data 
are used. It is therefore necessary to smooth the data, e.g., according to Eq. (2.12), 
before inserting them into Eq. (2.36), and to avoid a region close to the fleld boundary 
which can, however, be chosen small. If these precautions are taken, the flnite-fleld 
reconstruction is promising, especially since it can be combined with the non-linear 
reconstruction method. 

d. Further substantial improvement of the quality of the mass reconstruction is achieved 
by the synthetic fleld extension. The reason for this is simply that it reduces R in 
the effective fleld, and thus raises the lower bound on k, provided by the flnite-fleld 
reconstruction. As can be read off from Fig. 5, the uncertainty of the flnite-fleld mass 
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reconstruction relative to the average after field extension is ~ 0.1, and even close to 
the boundary the average reconstructed mass fraction is ~ 0.8. 
e. The mass reconstruction within a set of nested apertures, which also underestimates 
the true mass, performs equally well as the non-linear finite-field reconstruction after 
field extension in terms of its uncertainty, but it provides a less tight lower bound 
than the latter. 

On the whole, two main conclusions can be drawn. First, cluster reconstruction 
techniques which require a convolution integral to be performed over the entire real plane 
do not yield reliable mass reconstructions, because they introduce boundary effects which 
are hard to control. Overestimates of the mass are more likely and more substantial than 
underestimates, especially if the lens effect of the cluster under consideration is weak. 
Second, modified convolution kernels which apply to a finite field of integration yield lower 
bounds to the total mass. The tightness of these bounds depend on the size of the field. 
For the side length of 5' of the fields numerically simulated here, the reconstructed mass 
fraction /(r) drops from ~ 0.8 at the cluster center to ~ 0.3 at a distance of 2.'5, and 
this can be improved by synthetic field extension to /(r) ~ 1.0 at r = and /(r) ~ 0.8 
at r = 2.'5. Based on these results, it appears probable that the cluster reconstruction 
techniques based on the weak lens effect can achieve an accuracy of the reconstructed 
mass on the order of 10% . . . 20%, if they are applied to fields with side lengths of 10' or 
larger. 

Some cautionary remarks must finally be added. First, seeing deteriorates the intrin- 
sic signal, and guiding errors, telesope aberrations and similar effects can cause spurious 
signals. However, the exceptional quality of data taken with the refurbished HST (as an 
example, see plate Lll of Dressier et al. 1994), combined with very elaborate techniques 
of data analysis (cf. Kaiser et al. 1994b), give rise to the expectation that these difficulties 
can be overcome. Second, the redshift distribution of the arclet sources is yet unknown, 
and this introduces an uncertainty which can become substantial for high-redshift clusters. 
Third, all methods described here suffer from the scaling invariance (3.5) of the conver- 
gence, which leaves the observable signal unchanged. An additional sheet of constant 
surface density would therefore never be detectable with methods based on shear alone. 
This degeneracy can be removed by either requiring that «: — )• outward, or observation- 
ally by methods which are sensitive to the magnification of images. Fourth, clusters are 
embedded into large-scale structures, which by themselves contribute shear and magnifi- 
cation on a level of 10% . . . 20%. Possible effects of these additional large-scale lenses will 
be the subject of a further study. 

References 

Bartelmann, M., 1994, A&A, in press 

Bartelmann, M., Weiss., A., 1994, A&A, 287, 1 

Bartelmann, M., Steinmetz, M., Weiss, A., 1994, A&A, in press 

Dressier, A., Oemler, A., Sparks, W.B., Lucas, R.A., 1994, ApJ, 435, L23 

Kaiser, N., 1992, ApJ, 388, 272 

Kaiser, N., 1994, preprint CITA-94-39 or astro-ph/9408092 
Kaiser, N., Squires, G., 1993, ApJ, 404, 441 



24 



Cluster mass estimates from weak lensmg 



Kaiser, N., Squires, G., Fahlman, G.G., Woods, D., 1994a, preprint CITA-94-40 
Kaiser, N., Squires, G., Broadhurst, T., 1994b, preprint astro-ph/9411005 

Kaiser, N., Squires, G., Fahlman, G.G., Woods, D., Broadhurst, T., 1994c, preprint astro-ph/9411029 

Kochanek, C.S., 1990, MNRAS, 247, 135 

Miralda-Escude, J., 1991, ApJ, 370, 1 

Narayan, R., Bartelmann, M., 1994, preprint 

Schneider, P., 1994, preprint 

Schneider, P., Seitz, C, 1994, preprint 

Schneider, P., Ehlers, J., Falco, E.E., 1992, Gravitational Lenses (Heidelberg: Springer) 

Seitz, C, Schneider, P., 1994, preprint 

Tyson, J.A., 1988, AJ, 96, 1 

Tyson, J.A., Seitzer, P., 1988, ApJ, 335, 552 

Tyson, J.A., Valdes, F., Wenk, R.A., 1990, ApJ, 349, LI 

Webster, R.L., 1985, MNRAS, 213, 871 



